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I. INTRODUCTION 



A longstanding puzzle of general relativity is whether or not it is possible to forbid the 
presence of exotic phenomena, such as the construction of closed timelike curves ^ and 
superluminal travel |2]. A variety of energy conditions exist that would prevent these pos- 
sibilities. The weak energy condition (WEC) and null energy condition (NEC) respectively 
require that T^i^V^V^ > for any timelike and null vector V^, but these conditions have 
well-known violations even in free field theory in Minkowski space. 

The conditions may be weakened by requiring that they hold only when averaged over a 
complete geodesic, and those averaged energy conditions are still sufficient to rule out exotic 
phenomena. The averaged weak energy condition (AWEC) is violated in the standard elec- 
tromagnetic Casimir energy between parallel, perfectly conducting planes, and it appears 
that the violation occurs even with more realistic materials j^. The system does not vio- 
late the averaged null energy condition (ANEC), however, because geodesies in directions 
which violate NEC all pass through the conducting planes. Assuming any realistic material 
conductor in these places, the positive stress-energy there more than cancels the negative 
energy between the planes. In j3| it was shown that for a simple model of scalar fields, a 
static domain wall violates WEC, NEC and AWEC, even when the domain wall energy is 
included, but again ANEC is obeyed, for the same reason as above. 

ANEC is violated in curved space-time [1, H, 0,0]' but it is not clear that such violations 
can be used to produce exotic phenomenon. They are typically overwhelmed by ordinary 
gravitational energy densities, for objects that are not at the Planck scale. No violation of 
ANEC is known in a consistent quantum field theory in fiat spacetime. 

One can try to avoid the problem of passing through the material in the examples above 
by considering a geodesic passing outside a ball, where one can safely ignore the contribution 
to the energy density from the material of the ball. For the case of minimally coupled scalar 
fields, outside a perfectly refiecting spherical shell, one finds the same results as in the 
domain wall case: WEC and NEC are violated, but ANEC is obeyed 

In this paper, we analyze Casimir energy densities in the electromagnetic field outside a 
ball with fixed dielectric constant and permeability. The results also extend to the case of 
a perfectly conducting sphere.^ In contrast to the scalar field case, all energy conditions, 
including NEC (and therefore ANEC) are obeyed everywhere outside the ball. Thus there 
is still no counterexample to the conjecture that all quantum field theories in fiat spacetime 
obey ANEC. 

The generality of our result — depending only on the requirement that the index of 
refraction n is bigger than 1, which is necessary to maintain causality in classical electro- 
magnetism without absorption — suggests that it might be possible to extend it to other 
cases in QED or other conformal theories.^. 

A number of related calculations have been carried out in other contexts. The total energy 



^ In the case of finite, nonzero conductivity, we have absorption, couphng the incoming waves to modes 
within the material. The effects of this interaction might depend in detail on the internal dynamics of the 
material, so we do not consider it here. . . 

^ In the case of a conformal scalar field, using the techniques of |j9| one finds that outside a perfectly 
reflecting spherical boundary, NEC is also always obeyed. Unlike the electromagnetic case, however, this 
result does not hold in each channel individually; for the £ ^ partial wave contributions the NEC is not 
always obeyed, but when one includes all £ it is obeyed everywhere. 



2 



of a dielectric ball has been studied, in particular with applications to the controversial 
question of its relevance to sonoluminescence 0, 11 1. However, it remains difficult to 



formulate precisely because the energy is sensitive to the details of the dynamics within 
the material. These potential divergences were studied in detail in l3|; see also reviews 



|13L 114 |15[, and references therein. The case of the Casimir energy density outside a 
conducting sphere is a standard calculation, originating with the work of Boyer and 
reviewed in |il7|. Milton 3] calculates the total energy including the interior and the net 
force on a dielectric sphere, and in 19^ 2^ uses a comparison of the the stress density inside 
and outside the sphere to find a Casimir contribution to the total stress on the sphere. 

In the next section, we give the energy density and the radial and tangential pressures 
in terms of the electromagnetic field components. In Section III we decompose the fields 
in vector spherical harmonics, and in Section IV we derive the boundary conditions for the 
ball with constant dielectric constant and permeability. In section V we quantize our fields, 
and in Section VI we compute the energy density and pressures in the vacuum outside the 
ball. In Section VII we prove that all the energy conditions are obeyed, and we conclude in 
Section VIII. 



II. THE ELECTROMAGNETIC STRESS-ENERGY TENSOR 

For the electromagnetic field outside a ball of radius R, with constant dielectric constant 
and permeability, we will calculate the vacuum expectation value of the energy density, 
radial and tangential pressures. We will use Heaviside-Lorentz units throughout. 

The symmetric electromagnetic stress-energy tensor is 

j.,u ^ g^^F^^F^'' + -^g^^F^^F'^^ (1) 
The 00 component is the electromagnetic energy density, 

TOO = 1(^2^52) (2) 

and the stresses are 

= -{E'E^ + B'B^ - ]^5'^{E'^ + 5^)) (3) 

where i,j are Cartesian coordinates in three space. We can rewrite the diagonal components 
in terms of a spherical coordinate system, which we will decompose into radial and tangential 
components, denoted by subscripts r and t respectively. Let 

Fr = \{El + Bl) (4) 



Ft = \{E^, + B^,) . (5) 



Then the energy density and the radial and tangential pressures are respectively given by 

p = TOO = + (g^) 
Pr = = Ft-Fr (6b) 
Pt = T't"*' = Fr . (6c) 
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III. CLASSICAL FIELDS 



We will work with the electromagnetic potential A'^{'x,t). In empty space we can work 
in the radiation gauge, where d^j^A^ = and ^4° = 0, so that 

and 

B = VxA. (8) 

We separate the two polarization channels, so that we can write the vector potential as a 
sum of classical modes, 

v-^ f°° dk 

A(x, t) = 22 -7k= {Kim(^~''"'(^tem + Kim^~'''''\i^ + complcx conjugate) (9) 

where uj = k outside the ball since the field is massless, the sum over £ and m gives the 
partial wave expansion in the 3 spatial dimensions, and we define 

A+„(r,r]) = X,^(l])^+(r) (10) 
K,mir.^) = X (X,„(r])^,-,(r)) (11) 

in terms of the vector spherical harmonics 

The modes A^^^ are the magnetic multipoles, also known as the transverse electric (TE) 
modes, and A^^^ are the electric multipoles, or transverse magnetic (TM) modes. The 
radial functions ^'/^^(r) and ipj^ii^r) satisfy the same equation of motion as a scalar field, 

+ ^keir) = k^i'Ur) ■ (13) 



dr"^ r dr ^2 ^ 

For any function gkeir) obeying this equation, we have the following useful identities: 

f •X^^c/fc^(r) = (14) 
f • (V X (X,„^/fc,(r))) = -v/^(£+l)^7fc,(r)F,„ (15) 



r 

V X (V X (Xerr^gkeir))) = -V^ (X^„c/H(r)) = k^Xe„^gH{r) (16) 

1 d 

f X (V X (Ki^gkeir))) = — — (r^H(r))X<.^ . (17) 

r or 

The E and B fields, which are related to A via Eqs. ((Tj) and (jS)) can then be expressed 
771^ (E^,„e-'"*a+^ + E^.^e^'^^'a^,^ + complex conjugate) (18) 
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and likewise for B, with 



^tim = ^^^em^tiir) (19) 
B+„ = V X (X,„^+ (r)) (20) 

^kim = ^^^V X (X,„^,- (r)) (21) 



V X (^V X (X,^^- (r))^ = A;X,^^- (r) . (22) 



Using the identities above, we can now decompose the modes separately into radial and 
tangential components, 

^Mm,r = (23) 

K^^^t = i^^imi^Ur) (24) 

Btim,r = \4W^)yin.i^l,{r) (25) 
Kim, = XX,^ (26) 



and 



E^m^r = -^VW+^)Yem^;eir) (27) 

Efcw = ^^K".W)rxX,„ (28) 

BMm,r = (29) 

^kim,t = kXemip-^ir) . (30) 

IV. BOUNDARY CONDITIONS FOR THE UNIFORM BALL 

The general solution to the radial equation of motion, Eq. (fT^ . can be written 

i^Ur) = ^k''"^'^hf\z) + hf\z)] (31) 



/27r 

where z = kr, 5f{k) is the scattering phase shift, and 

,(1>2)/ N _ . A^7t(1,2) 



^y"^(^) = i/^^s;y^) (32) 

are the spherical Hankel functions. 

We now consider a uniform ball of radius i?, dielectric constant e, permeability /i, and 
zero conductivity. Nonzero and finite conductivity leads to absorption, in which case the 
Casimir energy may depend in detail on the properties of the conducting material, so we 
will not consider that possibility. The case of perfect conductivity, however, can be obtained 
as the limit where e — cx). 

Within the ball, we have the same decomposition of the fields, except there we have 
k = nuj, where n = y/ql is the index of refraction. We will consider only cases with index 
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of refraction n > 1, since a frequency- in dependent index of refraction below 1 would violate 
causality.^ We can also scale the wavefunction inside the ball by an overall constant in 
each channel. At the surface, eEr, Ej, Br, and -Bf are continuous, giving two independent 
conditions for the two unknown quantities: the phase shift and the normalization constant 



of the interior wavefunction. Solving for the TE phase shift, we have |21 

g2i5+(fc) ^ aj/{nxyhf\x) -Ji{nx)hf^'{x) ^^^^ 
aji {nx)h^P {x) — Ji{nx)h^P ' [x] 

where prime denotes differentiation with respect to the function's argument, and we have 
defined a = \/e/iJ,, x = kR, and the Riccati-Bessel and Riccati-Hankel functions 

/ TV Z 

n{z) = zjfXz) = J — Je+i/2{z) (34) 
J^^(z)^zl,r(z) = MH^,(z). (35) 



Defining 



we can write 



and thus 



six) = W\ (36) 



as{nx)hf\x) — hf^'{x) 



(37) 



^2iS+{k) _ ^ ^ n as{nx)je{x) -j/{x) ^^^^ 
as{nx)hf\x) — hf^'{x) 
The TM mode is the same with and S^{k) S^{k). 

V. QUANTIZATION 

To quantize the electromagnetic field we simply declare that the coefficients af^^ in Eq. © 
are operators and take the Hermitian conjugate of those operators in the complex conjugate 
term. For the operators to have the usual relations in which all commutators vanish except 

[akim, al'f'm'] = S{k- k')5ie5mm' ■ (39) 
We require the radial wave functions be normalized by 

fOO 

r'^|J^eiry^P^,,ir)dr = 6ik-k') (40) 







^ A real material can have n < 1 at some frequencies, but such materials necessarily have absorption as 
well. 
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which is satisfied by Eq. ()31|). The vector and scalar spherical harmonics are normalized by 
the completeness relations 



iX^ml^ — \yem\ 



,2 2£+l 

4:71 



(41) 



The quantum field A(x, t) then obeys the conventional equal-time commutation relations in 
radiation gauge, 



[Ai{^,t),Ej{y,t)] 



-I ( - ^ 



- y) . 



(42) 



Since we are interested in computing the vacuum expectation values of the components of 
the electromagnetic field, from this point on, notations such as and will denote the 
quantum mechanical vacuum expectation values of those quantities. 
For the TE mode we find 



El, 



Blr 





oo 

E 

1=1 

oo 

E 

£=1 

oo 

E 

^=1 



2£+ 1 



57r 



dkuj k^^(r)| 



2^ + l^{^ + l) 

2£+l dk 



Svr 







ld_ 

r dr 



(43) 
(44) 

(45) 

(46) 



The renormalized contribution is found by subtracting the free wave contribution, 
yielding 



Blr 


= 


El 


oo 

= E 






Blr 


oo 

= E 






Bl 


oo 

= E 

£=1 



2e + i 



Svr 

2i + l £{£+!] 
Stt 

2e + i r dk 



67r 



UJ 



dk 




2 


Jo ^ 












2 


\ r or 





-(0), 



(47) 
(48) 

(49) 
(50) 



Using the general form of the wavefunction, Eq. (jHH), and subtracting the free wavefunc- 
tions given by Eq. (j31|) with 6i = 0, we find 



27r 



-l]h'''\z?+{e 



-2i<5? 



-i]hf\zr 



and 



2Tir' 



(51) 



(52) 
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where prime denotes differentiation witli respect to z = kr, and we write Se instead of S£{k) 
for simplicity. The Appendix of showed that the second term of Eq. ()5H) is just the first 
term with the replacement k — >■ —k + ie and likewise for Eq. Thus we can drop the 

second terms in both equations and extend the range of integration over k to — oo, with the 
understanding that k is to be taken above any branch cut on the negative real axis, 





= 


Elt 


oo 

= E 




£=1 




oo 

= E 






Bit 


oo 

= E 

i=l 



(53a) 

dkook^ (e^'^^ -l)hf\zf (53b) 



2i+l 
167r2 

^£(£ + 1) r (e-^^ - l) hf\zr (53c) 
2i+l r dk^2f2.si _Aj^ii),^^^2 (53^) 



The calculation for the TM modes is exactly analogous 



=1 



S!, = (54c) 
B-t = E ^ (^'^'^' - l) ^J'^(^)' (54d) 



VI. ENERGY DENSITY AND PRESSURE 

Substituting Eqs. (jSSl) and ^ into Eqs. dH) and © we find 



1 roo 7,2 

= 1)^(^ + 1) / dkn-hf\zr (55) 

£=1 °° 

^* = ^^^2 + 1) / + -^i'^'(^)' (56) 

O/TT r J-oo L J 



where 



Together, these expressions reproduce the energy density found in |17| . 

Following the methods used in |^ and we now convert this expression to a contour 
integral, which we close in the upper half plane. Since our interaction has finite range and no 
bound states, the S*- matrix in Eq. (jHTf) is analytic everywhere in the upper half-plane |2l|. 
Thus the only contribution to the integral comes from the branch cut along the imaginary 
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k axis coming from lu = y/W+Je . To the right uj = \/k^ = k, but to the left u = —k, so 
with k = tK, we obtain 

-| />oo 

(58) 
(59) 



— ' poo 

— -,5^(2£+i)/ d.f.n{hicr-k/icr} 

i=i 



where ( = nr, 



U0 = Ch{0 = \lY^i+i/2{0 (60) 



is the modified Riccati-Bessel function of the third kind, and 
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From Eq. (jHHjl . we get 



= + ^ (62) 

afinx)h{x) - k/ix) 

where x = ^-R? 

fix) = s{ix) = 1^ (63) 



and 



MO = = J^h+i/2{C) (64) 



is the modified Riccati-Bessel function of the first kind. 
VII. POSITIVITY 

We would like to know the sign of and Ft. First consider F^. Everything on the right 
hand side of Eq. (|3S|) is manifestly positive, except for and the overall negative sign, so 
Fj. has the opposite sign from T),. 

Now consider Ft. The properties of modified Riccati-Bessel functions are discussed in 
Appendix ^ From Eq. ()A12|) we find 

fc/(C) < -HQ (65) 

Since both sides are negative, 

kaCf > UC? (66) 

so the term in braces in Eq. (j59|) is negative, and Ft also has the opposite sign from T^. 
Let us now determine the sign of T^. We can write in the form 

aa — h a^^a — b 2ac + 2bd — {a + a~^) {ad + be) , . 

ac — d a~^c — d + d'^ — (a + a~^) cd 
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with 



« = f{nx)%{x) (68) 
b = t/ix) (69) 

c = fMHx) (70) 
d = k/ix) (71) 

The functions % and ki are positive, and % is increasing, so / is positive and a, 6, and 
c are positive. On the other hand, ke is a decreasing function, so is negative. Thus each 
term in the denominator of Eq. (j67|) contributes positively. 

We can write the numerator of Eq. (jU7j) as 

2(a-b){c- d) - {a + a'^ - 2) (ad + be) (72) 

and we have 

ad + bc = f{nx) %{x)Ki.x) + h{.X%{.x) (73) 

By multiplying Eq. ()A14|) by the positive quantity %{x)^f.{x) we find that the term in 
brackets in Eq. (fTSj) is positive. Since f{nx) > 0, and {a + a^^ — 2) = (a — lY/a > 0, the 
second term in Eq. (ff^ contributes negatively. 
Now c — c? is manifestly positive, and 

a-b=%ix)[finx)-fix)] (74) 

Since n > 1 by assumption, Eq. ()A6|) tells us that the term in brackets is negative, so the 
first term in Eq. (f7^ is negative. Thus < 0, and both Fr and Ft are positive. 

From Eqs. (jH)), we see that p > pt > and p > \pr\- Thus the weak, null, dominant, 
and strong energy conditions are satisfied at each point, and therefore the corresponding 
averaged energy conditions are satisfied on every geodesic that lies outside the ball. 



VIII. DISCUSSION 

We have demonstrated, without resorting to numerical calculations, that the electromag- 
netic Casimir energy density is positive outside a ball with constant dielectric constant and 
permeability, and that all the standard energy conditions are obeyed. We needed to assume 
only that the index of refraction is greater than 1, which is necessary to maintain causality. 
We did not consider the cases of absorption and dispersion. We expect, however, that our 
results would hold in such cases, and furthermore, we expect them to remain valid in the 
case of material properties that vary with radius. 

The electromagnetic situation discussed here gives a simpler result than the case of a 
scalar field. For the minimally coupled scalar field |9], the pointwise energy conditions are 
not obeyed, although the average null energy condition nevertheless holds for every geodesic 
outside the ball. For the conformally coupled case the null energy condition does hold at each 
point, but only when all angular momentum modes are considered. In the electromagnetic 
case, the energy conditions hold at each point separately for each k and /. 
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APPENDIX A: RICCATI-BESSEL FUNCTION INEQUALITIES 

We let Ti and k£ be modified Riccati-Bessel functions as in Eqs. (jHH) and (plj) . and define 
their logarithmic derivatives, 

/W=S^ (Al) 

as in Eq. (jU^ and 

9{x) = (A2) 
Now the modified Riccati-Bessel functions satisfy the differential equation 

x^r; - (x^ + i{i + 1)) j-^ = (A3) 

with JF = Tor k, so the derivatives satisfy a generalized Riccati equation (j^ . 



for h = f or g. 

Differentiating Eq. ()A4|) gives 



h"(x) = -2^^^^ - 2h(x)h'(x) (A5) 

X'^ 

from which we conclude that if h'{x) = for some x then h"{x) < there. Thus h'{x) 
cannot increase through 0, and so if h'{x) is negative for some x, then it is negative for all 
larger x, while if it is positive for some x, then it must be positive for all smaller x. 
Now for a; — s> 0, f{x) ^ {I + and so for small x, f'{x) is negative, so 

fix) < (A6) 

for all X. As X — > oo, /(x) — * 1, so 



for all X. 

For X — >• oo, 



/(x) > 1 (A7) 



so 



k,{x) -e-^ [I + + 0{x-') ) (A8) 

%{x) - -U^) - \e-^ (^^^ + 0{x'-') ) (A9) 



11 



therefore 

( ^ 1 ^( 

nl T I — > — I — — 

2a;2 



g{x)^-l-^-^^ + Oix-') (AlO) 



Since g'{x) > for large x, we have 

g'{x) > (All) 

for all X, and thus 

g{x) < -1 (A12) 

for all X. 

Now by subtracting Eq. ()A4jl for / from the same equation for g, we find 

g>(x)-f{x) = f{xf-g{xf (A13) 



From Eqs. (jA6jl and (jAllll this quantity is positive, so /(x)^ > g^x)"^. Since (^(x) < while 
f{x) > 0, we find 

f{x)+gix)>0 (A14) 
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